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Abstract 



Q^ . We present in this paper the general framework of a method which permits to restore 

^^ \ the rotational and particle number symmetries of wave functions obtained in Skyrme 

HF+BCS calculations. This restoration is nothing but a projection of mean-field in- 
trinsic wave functions onto good particle number and good angular momentum. The 
O \ method allows also to mix projected wave functions. Such a configuration mixing 

is discussed for sets of HF-I-BCS intrinsic states generated in constrained calcula- 
tions with suitable collective variables. This procedure gives collective states which 
are eigenstates of the particle number and the angular momentum operators and 
^ ' between which transition probabilities are calculated. An application to ^^Mg is pre- 

j^ ■ sented, with mean-field wave functions generated by axial quadrupole constraints. 

Theoretical spectra and transition probabilities are compared to the experiment. 

PACS numbers: 21.10.Ky 21.30.-n 21.60.Jz 27.30-Ft 



1 Introduction 



The cranking method is widely used in nuclear spectroscopy to describe high 
spin states. Applications based on effective nuclear interactions have been par- 
ticularly successful in the description of super-deformed rotational bands in 
several regions of the mass table [1-5]. In the cranking method, a rotational 
band is generated by the rotation of a deformed intrinsic state. Since cranking 
states are not eigenstates of angular momentum, it is not straightforward to 
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determine transition rates in nuclei which are not very well deformed. To over- 
come this difficulty, approximations [6-8] have been developed for transitions 
within a band. However, they are only valid when the structure of the nuclear 
states are not affected by rotation, i.e. for rigid nuclei. 

Another limitation of the cranking model occurs in nuclei soft with respect 
to the variation of a collective variable. In this case, one expects that the 
interference of the zero-point vibrational mode with the rotational motion 
will lead to variations in the nuclear structure along the yrast line. 

There have been several attempts in the past to introduce phenomenological 
corrections taking into account rotational and vibrational corrections to the 
mean-field energy (see for instance [9-11]). It is the aim of the present paper to 
introduce a method in which rotations and vibrations are taken into account 
simultaneously in a more general and consistent way. 

The starting point of our approach is a set of many-body wave functions 
generated by constrained Skyrme HF-I-BCS calculations. The discretization of 
these wave functions on a 3-dimensional Cartesian mesh enables to describe 
very general shapes of the nuclear density and to easily write the effect of a 
spatial rotation on the mean-field wave functions. This property permits to 
restore symmetries with respect to angular momentum and to the proton and 
neutron particle numbers in a systematic way. 

It is not the ffist time that symmetries are restored on mean-field wave func- 
tions. Several methods have been developed for schematic nuclear interac- 
tions and a limited number of active valence nucleons (for recent references, 
see [12,13]). Methods based on effective interactions have been limited to ei- 
ther angular momentum [14-16] or particle number projection with mixing of 
configurations by the generator coordinate method (GCM) [17]. The present 
work generalizes this last work by the inclusion of a projection on angular 
momentum. 

We present below the general framework of our method together with a test 
on a light nucleus for which extensive calculations can be performed. In the 
first part, we show how to implement the projection on J, A^ and Z simul- 
taneously and to determine the contractions needed to calculate the matrix 
elements between two different many-body HFB vacua. The specificities of 
mesh calculations require to take explicitly into account the fact that the 
single-particle bases are not complete and lead to formulae rather different 
from some similar works [18,13]. In the second part of this work, we present 
an application to the ^^Mg nucleus for which the properties of several effective 
interactions are tested. 



2 Angular momentum and particle number projections 

2. 1 Principle of the method 



The starting point of our metliod is a set of wave functions \^a) generated by 
mean-field calculations with a constraint on a collective coordinate a. Wave 
functions with good angular momentum and particle numbers are obtained 
by restorations of symmetry on |$a): 

\<^,JMa) = ^Y.9KpitKP'P''\^a) , (1) 

where A/" is a normalization factor, and the operators P are projectors. 

A configuration mixing on the collective variable a is then performed for each 
angular momentum: 

|vl/,JM) = ^/i^^|$,JMa) . (2) 

a 

The weight functions /^^ are found by requiring that the expectation value 
of the energy: 

^j^ ^ {^,JM\H\^,JM) 

{m,JM\m,JM) ' ^ ' 

is stationary with respect to an arbitrary variation 5f^^. This prescription 
leads to the discretized Hill- Wheeler equation [19]: 

Sr^f'uJM rpJMTJM\ rJM,k r, / a\ 

a 

in which the Hamiltonian kernel Ti"^^ and the overlap kernel X"'^ are defined 
as 

ni^, = {<^JMa\H\(^JMa'), I^^, = {^JMa\^JMa') . (5) 

Since the Hamiltonian is rotationally invariant (we will come back to the prob- 
lem of the density dependence of the interaction in section 2.8) and conserves 
the number of particles, one has to restore the symmetries on only one of the 
two wave functions entering in each matrix element like eq. (5). The kernels 
are obtained by integration on three Euler angles and two gauge angles of the 
matrix elements between rotated wave functions. 

Besides these kernels, we will calculate transition probabilities between differ- 
ent eigenstates of the Hill- Wheeler equation. This requires the calculation of 
the matrix elements of a tensor of order L, T/^, between projected states. 



Such a secular problem based on the configuration mixing defined by eq. 2 
amounts to a variation after projection in a many-body Hilbert space built on a 
limited set of states obtained for different values of the collective variables a's. 



2.2 Determination oj the contractions 



The matrix elements (eq. 5) are calculated in two steps. They are first deter- 
mined for a given rotation and for given gauge angles. Then they are integrated 
over the Euler and gauge angles. 

To compute the kernels, we make use of the Balian-Brezin theorem [20] which 
expresses them in terms of contractions of creation and/or annihilation single- 
particle operators. 

To simplify the presentation, the spatial rotation is applied on the left wave 
function, whereas the particle number rotation is applied on the right one. The 
wave functions after rotation are denoted by \L) and |-R) and the associated 
quasiparticle annihilation operators by /^, r^ respectively, with /i = 1, . . . , A^. 
The number N is the total number of HF states included in the mean-field 
calculations. 

The quasiparticle operators are obtained by the diagonalization of the Bogoli- 
ubov equations on a single-particle basis. When the equations are discretised 
on a 3D mesh, a convenient basis is provided by the eigenstates of the HF 
Hamiltonian. Since only a limited number of single-particle HF states are de- 
termined, the left and right bases are different. The particle annihilation and 
creation operators of the left and right particle bases will be denoted (a, a"*") 
and (6, b'^) respectively. They are related by a unitary transformation: 



a 



M = H Rt^,uK ■ (6) 



The quasiparticle operators /^ and r^ are expressed as a function of the left 
and right single-particle bases by Bogoliubov transformations [8]: 





WU ; (7) 



where the matrices Uf and V^ have dimensions N x N. The expressions for 
the right quasi-particles are similar. 



Using the same notations as in [21], we define the two vectors 



^hM . P^ ; 



in terms of which any two-body operator O can be expressed. 

The quasiparticle bases are related by a hnear transformation T. Taking into 
account that the matrix R hnking the HF single-particle bases does not mix 
creation and annihilation operators, one obtains: 

A = Tp, T = r ' , (9) 

\^ -E* D-^j 

where, 

D = {UP'^RUr + Vi^^RVr)-\ E = -{Ul'^RVr + V.^^RUD. (10) 



Balian and Brezin [20] have established the relation between T and the con- 
tractions (Ap), {pp) and (AA): 



/(ll) A I ED l\ 
(AA) = = , (11) 

oy V "^ 0^ 



/o 1 \ /o 1 , 

(pp) = = , (12) 

\^0 (r+r+)y 1^0 DE* 



/0(lr+)\ hD^\ 

<^^>io ojiooj ^ <"' 

In addition, they have derived a formula for the modulus of the overlap be- 
tween left and right many-body wave functions: 

{L\R) = ±{dei D-^f'^ (14) 

One obtains the contractions in the particle bases: 






(15) 



2.3 Symmetry restrictions 



The formulae derived in section 2.2 are general. In the present application, we 
impose several symmetry restrictions on the mean-field wave functions. These 
restrictions which we develop below do not limit the formalism, they are only 
intended to simplify the equations for our test case. They can be released, if 
need be, at the only cost of additional computational time. 

First, the total wave function is symmetric with respect to refiections across 
the X = 0, y = and z = planes and time- reversal invariant. With this sym- 
metry restriction to ellipsoidal configurations, the single-particle wave func- 
tions are described by 4-dimensional vectors [22], 






i?e$^(r,-,+,-p) 



(16) 



where the components describe, respectively, the real and imaginary spin up 
and down parts of the single-particle wave functions. The total parity is given 
by p. Each of the components has a well defined symmetry with respect to 
refiections across the three planes x = 0, y = and z = (see Ref. [21], 
appendix C). In a mesh calculation, the space is limited to a box beyond 
which all the wave functions are set to zero. With these symmetries, all the 
integrations can be limited to an octant of the box. 

A second restriction is that the intrinsic wave functions have been generated 
by constraints on the axial quadrupole moment only, the z-axis being the 
symmetry axis of the nucleus. To restore angular momentum requires then 
only to rotate the mean-field wave function by an angle (3 around the y-axis. 
The coordinates {xo,yo,Zo) become after rotation: 



xi = zq sin(/9) + xq cos(/9) 

1/1=1/0 , 

Zi = Zq COs(/5) — Xq sin(/5) 



(17) 



and the four components of a wave function rotate according to: 



Pi,^(ri) = w^i,A«(ro) cos(/9/2) - ^s.^l^o) sin(/5/2) 
P2,/.(ri) = W2,^(ro) cos(/3/2) - W4,M(ro) sm(/3/2) 
P3,M(ri) = W3,^(ro) cos(/5/2) + wi,^(ro) sin(/3/2) 
P4,^(ri) = u'4,^(ro) cos(/5/2) + W2,^,{ro) sin(/5/2) . (18) 

These equations show that while the total parity p is preserved by a rotation, 
the specific symmetries of the four components with respect to the a; = and 
z = planes are lost. In the same way, signature as defined in ref [22] is not 
conserved. An immediate consequence of this loss of symmetry is that matrix 
elements between rotated and non rotated wave functions must be calculated 
in half of a box. The real and imaginary parts of the rotated wave function 
have different symmetries with respect to the y = plane and the overlap 
between rotated and non rotated individual wave functions are real. 

Furthermore, the points {xi,yi,zi) do not coincide after rotation with mesh 
points. The values of the rotated wave functions on the mesh points are reeval- 
uated using the analytical forms of functions defined on a mesh given in ref [23]. 

The relation between rotated (left) and non rotated (right) operators is given 
by: 

<l^i = '^{Rfl,yK + Rfl,ubT7} , ajC = ^{Rjl^yK + Rji^yhy] , (19) 

where the state \jl) denotes the time-reversed partner of state |/i) and the 
matrix R is given by the overlap between left and right states. Using the 
properties of the time reversal operator, one can show the following relations 
between the overlaps: 

The matrix R does not separate into blocks corresponding to the different 
combinations of parity and signature as in our previous study [17]. 



2.4 Restriction to BCS transformations 



One can still simplify the expressions for contractions derived in sub-section 2.2 
for the general Bogoliubov transformations to BCS transformations. This lim- 
itation should not be too bad for even nuclei close to the stability line, as 
^"^Mg, for which pairing correlations do not couple significantly bound and 
continuum states. 

The state \R) is obtained from the original intrinsic state by multiplying by a 
phase e^'^'^ the occupation numbers in the original BCS (or HFB) transforma- 



tion 



\R) ^ |i?(0)) 



^r/i"/i '^T^(^ "ji 



UrtJ>p + t^r/.e*2<P6+ 



with Ura, Vra real and positive, and Ur-p = u 



''rjj.j '^r/i 



Vra = -n 



rfi' 



In the BCS case the matrices f/; , VJ , U^ , V^ take the form 



U^^Ur, V^^^-aVr, Ui 



U^ = Ur ; Vf = -aVre-'^'^ ; f/, = | ' \ ; Vr 




V, 




Vr 
Vr 



(21) 



(22) 



The matrices ui, vi, Ur, Vr, are real, diagonal and of dimension N/2 x N/2. 
The matrix a is defined as 



a 



-1 

1 



(23) 



Thanks to the symmetry properties given by equation (20), the matrices D 
and E introduced in subsection 2.2 have the same block structure as matrix 
R. These same symmetries allow the relation a^Ra = R. We can thus write 
the contractions, as well as matrices D, E, in a simpler way. We have: 



(a+b) = ViD'^Vre^''^ ; 
(a+a+) = -a{ViUi - ViEDVi) ; 



(bb) = a{UrVr + VrDE*Vre'"f)e 



2i<j)\ 2i<f> . 



D = ViRVre^'"^ + UiRUr ; 
E = aE] 

E = UlRVrC^'^ - VlRUr . 



2. 5 Elimination of non occupied states 



In the derivation shown in the above section, it is assumed that left and right 
bases are either complete or truncated in such a way that they span the same 
space. As discussed in ref [21], this property is not valid when orbitals are 
discretized on a three dimensional mesh. This problem can be solved by taking 
into account that the missing part of the expansion of the left states on the 
right basis includes empty states that do not affect the structure of the nucleus. 



These states are defined by the condition f ^ = and contribute neither to the 
overlap nor to the contractions. Since the structure of the matrices is different 
than in ref [21], we give now the formulae corresponding to the present case. 



Let us introduce: 



Vi=\ \ ; U,= \ ' \ ■ (24) 

0/ I 1 



and split the unitary matrix R in the form 



,7^ J, 
R=\ I . (25) 



With these notations, the contractions can be written as: 



(a+b) = I ' I ; (26) 

0, 



{hh) = a\ r . . I . ^27) 

0, 



(a^a^) = —a \ | . (2J 

0^ 

The matrices V, Sr and Si are defined by: 



V=(Ui{n+)-'Ur + VinVre^"^) ' ; (29) 



Sr = t/^(7^+)"V, - Vi'JZUre^'''' ; (30) 



Si = Ui'}ZVre^'^-Vi{n^)-'Ur . (31) 



2. 6 Determination of the overlap 



2. 6. 1 Modulus 



Equation (14) gives the overlaps up to a sign. Several methods have been 
developed to calculate the global phase of the overlap. We show in the next 
section which procedure we have implemented. In this section, we give the 
explicit formulae for the modulus of the overlap. We have: 



D-^ = UiRUr + ViRVre^''^ , 
so that: 



(32) 



{L\R) = ±{det{D-')) 



-Uxl/2 



± 



/ /_ _ _ _ _ \ \ 1/2 



det 



V 



V 



rUr u 

It is easy, although tedious, to demonstrate that the matrix D^ 



(33) 



D 



V: 



vui{n+)-^T+ 



\ 



(i+{Tz+y^UrV- u-^ + i+{Tz+)-^UrVUi{n+)-^T+ J 



(34) 



where V is given by eq. 29, can be expressed as the product of four matri- 
ces [24]: 






D 



u, o\ ln+ T+\ /«-! o\ n -Pr 

1 j \ J+ W+ j \ 1 j \ 1 



with: 



a, = Ui + VilZVrU, ^n+e^' 






(3r = ViUVrU, ^T+e^'^ ■ 



so that 



(35) 



(36) 






detL' = det(f/^ ) det(i?+) det(a;^) 
((7^+)-^)l det(P) , 



det(f//a;^) 



det 



(37) 



10 



where we have taken into account that the determinant of the unitary matrix 
R equals 1. We finally obtain: 



<^|''>(«-±(SS)"^ ■ P^ 



2.6.2 Phase 

For any space rotation angle f3, the overlap between two mean-field wave 
functions is positive for 0„ = and 0p = 0. We need to follow the evolution of 
the phase of the overlap along the integration paths on 0„ and (pp. One has: 



(L|i?)(0) = v'(Re[(L|i?)(0)])2 + (Im[(L|i?)(0)])V^W . (39) 

We can rewrite the global phase Q as: 

n(0) = Im \n{{L\R) (0)) = $(0) + n(0)7r , (40) 

where $, which can be defined as: 

$(0) = arctan(Im[(L|i?)(0)]/Re[(L|i?)(0)]) , (41) 

is limited to the interval [— 7r/2 , 7r/2]. The integer number n remains to be 
determined. The total phase and its derivative are given by: 



It can also be calculated directly from the matrix element of the particle 
number operator: 



{L\N\R)i<P) _ 

= 2 Y. {R,AKh')<^ + R,-A<hi')^} ^ (43) 

where we have used the symmetry properties discussed in the previous section 
and eq. (26) for the contractions (a+b). 
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Using this formula, we determine the phase of the overlap with a method 
similar to the one developed in ref [25,13]. Starting from = 0, for which 
n(0)=O, the value of the phase at neighboring angles is determined by: 



n(0 + 50)^^(0) + ^ 



A A 

— fi(0) + — fi(0 + 50) 



(44) 



where we have taken into account that the interval of integration has a length 
n and is divided into L segments of equal dimension 50 = tt/L. From the 
comparison between f2(0 + 6(f)) calculated from (44) and $(0 + 6(f)) obtained 
from (41), one determines n(0), provided the mesh is dense enough [6(f> small 
enough) . 



2.7 Calculation of multipole moments and transition probabilities 



The determination of transition probabilities requires the calculation of the 
matrix element of a tensor of order L, T/^, between eigenstates of the angular 
momentum operator. 

An eigenstate of the angular momentum operator, with eigenvalue J, is ob- 
tained by projecting the mean-field wave function |$): 



1$, JM) = ^K9KPi,K\^) ^43^ 

where the projector is given by [8]: 



2 T -\- 1 r 

Pmk = -^ J dnD'i^j,{n)R{n) , (46) 

with Q = (q;,/3,7) the Euler angles and R{Q) = e*"•^^e*^■^J'e*'^■^^ 

Only the K=0 term of this expression is present with the symmetry restrictions 
that we have imposed. Then: 



{JM, $|f 0| J'm', $') = (^l^oM^L^M'ol^) (47) 

where we have used the properties {P'mk)^ = P^m ^"^^ P'mkP'm'k' = ^.jj'^km'Pk 
of the angular momentum projector operator [26]. 

We have 
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(JM,$|f°|/M',<l>') 



^ {2J + 1){2J' + I) 



87r2 



!d^d9.'DiU^)Dil^{9.'){^R{^)TlR{9.')\^') 



-,1/2 • 



(4J 



The axial symmetry reduces the integration interval to [0,7r/2]. Moreover, 
thanks to the transformation of the wave functions under rotation (see eq. 18), 
the matrix element ($|T°e*^'^'' |$') is real and is the same for rotations of angles 
j3 and —fi. In the case of axial symmetry, the final expression takes then the 
form: 



(J'ML0|JM)(J'0L0|J0)5 



^ ' ' ^' ' ^ 2V 2J+1 



MM' 



1/2 



/sin/5ci/5ci^o(/5)(*|e^^^H*) J smp'dp'diQ{p'){<^'\e^f^'Jy\<^') 



1/2 • 



(49) 



In the case of electric quadrupole transitions, the diagonal matrix element 
takes the form: 



{JM = 0, ^\Q2o\JM = 0, $) = (J020| J0)2 
(J+1)J 



Jsmf3df3dUf3)me^^'^Q2o\<^) 



/sin/3rf/3rf4(/3)($|e^'3^H'^) 
Jsmf3df3d-loW){^\e^f'yQ2o\<^) 



(2J + 3)(2J- 1) Jsm(3d(3dU(3){^\e'^^'\^) 



(50) 



The transition matrix elements between GCM states are obtained as the 
weighted sums of the contributions of the different basis states. 

We finally obtain for the reduced transition probability between the initial, 
Ij-th GCM collective state of spin Ij, to the final, l/-th GCM collective state 
of spin If-. 



B{EL, Z 



(h) 



rih)^ 



{IiOLO\IfOy 



JJdada'f^^f'^f^*{a)f^^t'^^\a')Jdcos(3dii{(3){^{a)\e'^^yMLo\^'ia')) + {I+l 



f^'f 



-L „ .1/7, 



Udada'f^^^'^i^*{a)0'^f\a')Sdcospd^^M{'^{o)\^'^M^'{ci')) 



'Jf > h 



It, Id 



(51) 



with Alio = r'^YLo{0,^ 
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2. 8 Density dependence of the effective interactions 



The density dependent term of the interaction must be generahzed to calculate 
non diagonal matrix elements. In the case of a density dependence equivalent 
to a 3-body interaction, the Hamiltonians kernel can be expressed in terms of 
the left right mixed density [21]: 

p(r) = ^(a+b)^,$*^(r, a)$,,.(r, a) (52) 

liver 

We have chosen the same dependence on the mixed density when there is 
no equivalence with a three-body interaction. The energy is then expressed 
as a functional of |-R) and \L) similar to the mean- field functional. One can 
show that the mixed density depends only on the relative angles between the 
principal axes of \R) and \L). Therefore, after integration on the Euler angles, 
the energy is real and does not depend on the orientation of the reference 
frame. One can thus restore symmetries either on the left or the right wave 
function. 



3 Application to ^^Mg 



The results shown in this section have been obtained using the HF+BCS wave 
functions generated with an axial quadrupole constraint. The Lipkin-Nogami 
prescription has been used to improve the treatment of pairing correlations. 
It has indeed been shown that this prescription permits to generate wave 
functions which are reasonable approximations of those obtained by a variation 
after projection on the good particle number [27]. In this way, the lack of a 
complete variation after projection should be partly compensated. The mean- 
field results that we will present below correspond to these HF-I-BCS+LN 
calculations. 

Since our aim is mainly to test the properties of our method, we have per- 
formed calculations with three Skyrme parameterizations which have given 
satisfactory results in the description of rotational bands in well deformed 
nuclei, namely SIII [28], SkM* [29], and SLy4 [30]. The pairing interaction 
is a zero range interaction similar to the ones used in previous studies of 
super-deformed bands [2], and of nuclei far from stability [31]. In calculations 
performed with SLy4, we have varied the strength of the density-dependent 
pairing force from G = l250MeV fm^ to G = 900MeV fm^. Most results did 
not depend significantly on the value of G, the most sensitive quantity being 
the excitation energy of the first 2"*" state. We have chosen to show here only 
results obtained with G = lOOOMeVfm^ which gives a 2"*" energy close to 
experiment. The same value of the pairing strength has been used for the two 
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other Skyrme parametrizations. It is clear that a precise adjustment of the 
pairing strength requires a study of a large range of isotopes, a work which is 
in progress. We have also tested the density independent zero-range pairing 
interaction that has been adjusted for a study of super-heavy nuclei [32]. 

In figure 1 are shown the mean-field energies as a function of prolate defor- 
mations with the three Skyrme interactions and, in the case of SLy4, of a 
surface and a volume pairing force (dashed curves). Since oblate deformations 
play a minor role in this nucleus, they have not been represented. The curves 
are similar, with a well deformed prolate minimum corresponding to a mass 
quadrupole moment of approximately 1 b. The curvature is slightly different 
in the four cases, the volume pairing leading to the deepest minimum. 

The energies obtained by projecting each of the mean-field wave functions on 
good particle number and angular momentum are also shown on figure 1. The 
abscissa of the projected energies correspond to the quadrupole moment of 
the intrinsic wave function. The spherical configuration is a pure O"*" state and 



20 



> 

0) 



LU 



10 
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Sly4, Vol. pai. 



(b) 



(c) 
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Sly4, G=1000 SIII,G=1000 SkM*, G=1000 











2 
qo(b) 







Fig. 1. Mean-field and projected energies obtained for ^^Mg as a function of the axial 
quadrupole moment qo. The pairing interaction is a zero range interaction without 
(a) and with (b-d) a density dependence. The Skyrme parameterizations are SLy4 
(a and b), SIII (c) and SkM* (d). The dashed hues correspond to the HF-^BCS-hLN 
energies. The energies obtained by projecting on angular momentum (from 0"*" to 
8"*") intrinsic wave functions curves are plotted in fuh hne as a function of the 
quadrupole moment of the intrinsic wave function. 
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contributes only to the 0^ projected curve. The energy gained by projection 
in this case is due to the difference between the Lipkin Nogami approximation 
of the energy gain due to projection on particle number, — A2AA^^, and the 
exact gain. It is of the order of 1 MeV in all cases, except in the calculation 
with the Sill force, where it is significantly larger. This may be due to that a 
lower value for the pairing energy is obtained in the Sill calculation and that 
the LN prescription is less accurate in the low pairing regime. 

In all cases, the projection on angular momentum increases the energy differ- 
ence between the spherical configuration and the minimum of the J=0+ curve 
by 3 to 4 MeV. The intrinsic wave function leading to the minimum of this 
curve has a quadrupole moment slightly larger than the one corresponding to 
the mean-field minimum. For higher angular momenta, the minima are shifted 
to higher quadrupole moments. The dependence of the curves on the nuclear 
interaction is rather weak for all J values other than 0; for J=0, the differences 
are mainly related to differences between the spherical configurations. 

On figure 2, the weights of the various angular momenta in the mean-field 
wave functions are plotted as a function of the axial quadrupole moment. 
These weights do not depend significantly on the nuclear interaction; results 
obtained with the SLy4 force and a surface pairing are only shown. 

The spherical configuration is already an eigenstate of the angular momentum 
operator and is affected only by the particle number projection. Its N=Z=12 
component has a weight of 0.18, the remaining part corresponding to differ- 
ent proton and neutron numbers. The J=0 curve displays a pronounced peak 
around the spherical point. It is not symmetric with respect to quadrupole 
moment, the weights on the prolate side being larger than on the oblate side. 
For this reason, the J=0 curve is rather far from a Gaussian. This shows 
that the Gaussian overlap approximation sometimes used to derive a collec- 
tive Shrodinger equation should be used with care. When the weight of a 
mean-field wave function component is lower than 0.01, the projection be- 
comes numerically difficult; such tiny components have been excluded from 
the configuration mixing calculations. As expected, the quadrupole moment 
correspondig to the largest weight moves to higher quadrupole moments as a 
function of angular momentum. For J=10'^, no maximum is obtained up to a 
quadrupole moment of 4.5 b, showing that our variational space is inadequate 
to describe accurately states with such a high spin. A dissymetry between 
the prolate and oblate deformations can also be noticed on figure 2: the total 
weights corresponding to N=Z=12 are twice larger for prolate configurations 
than for oblate ones. This confirms that ^^Mg is dominated by prolate config- 
urations. The different dispersions in particle numbers for prolate and oblate 
deformations is at the origin of this dissymetry (upper part of the figure). 

The variation of the energy as a function of prolate and oblate deformations is 
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Fig. 2. Lower part: weights (^j''^'Lol'^(9o)) as a function of qo and J calculated 
in case (b) of fig. 1. Upper part: particle number dispersions in the intrinsic wave 
functions as a function of go for neutrons (full line) and for protons (dashed line). 



17 



plotted on figure 3 for tlie SLy4 Skyrme parametrization and a surface pairing 
interaction. In addition to the prolate minimum already observed on figure 
1, the mean- field curve presents a shoulder at an oblate deformation around 
0.5 b. The full projection creates an oblate minimum at the position of that 
shoulder for J values ranging from 0^ to 6^. For greater values of J, the weights 
of the intrinsic wave functions for deformations below -2 b are very small. 
Consequently, the projected energy curves do not exhibit any oblate minima. 
However the J=0"'" to 6~^ minima are probably not stable against triaxial 
deformations, since a calculation including triaxial deformations indicates that 
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Fig. 3. Projected energies for ^^Mg as a function of the axial quadrupole moment 
(case (b) of fig. 1). The symbols are the same as in the previous figures. The first 
three energies obtained for each angular momentum in the configuration mixing 
calculation are represented by horizontal bars centered at the value of qo where the 
respective collective wave functions are maximum. 



the shoulder in the intrinsic curve is a maximum with respect to 7. 

For each value of the angular momentum, we have performed a configuration 
mixing calculation including quadrupole moments between -3.5 b and 4.5 b. 
This corresponds to intrinsic configurations excited by about 30MeV with 
respect to the prolate minimum. This configuration mixing is nothing but a 
variation after projection in a limited but hopefully relevant Hilbert space. 
The spectrum generated in this way (represented by bars) is plotted at the 
quadrupole moment corresponding to the largest component of the collective 




Fig. 4. Squared weights of intrinsic states in the Yrast levels (/ 
angular momentum representations. 



0) of the different 
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wave function, which are shown on figure 4 for the yrast states. The value 
of this quadrupole moment is very close to the minimum of the projected 
energy curve. Moreover, the energy of this minimum is slightly modified by the 
configuration mixing. The largest gain, ~800 keV, is obtained for the 0+ state, 
but is reduced at higher spins. Several excited states are found at low energy 
for each spin value. Except for the second 0^ and 10^, the wave functions of 
yrare states are peaked around the oblate secondary minimum. 



The spectra obtained with the four choices of interactions described for fig. 1 
are compared to the experimental spectrum [33] on figure 5. The energies cor- 
respond to the minima of the projected energy curves except in the column 
(b') where the configuration mixing spectrum is represented. For compari- 
son, we show also the results of a cranking calculation performed with the 
SLy4 interaction and a surface pairing with the same strength of 1250MeV 
that we used in our calculations of SD rotational bands. In this case, we use 
the HFB+Lipkin-Nogami method presented in our previous cranking calcula- 
tions [2]. 

Since the intrinsic ground state of ^'^Mg has a large prolate deformation, one 
expects the cranking approximation to be valid. The cranking spectrum is 
indeed in good agreement with the experimental data up to the 6~^ level. Ex- 
(Exp.) ^ (Crank.) ^ (a) ^ (b) ^ (b') ^ (c) ^ (d) 
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Fig. 5. Experimental positive parity levels of ^^Mg compared to a cranking cal- 
culation (SLy4 and density-dependent zero-range pairing interactions) and to the 
(a)-(d) choices of interactions described in fig.l. In the column (b') is plotted the 
configuration mixing spectrum. 
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perimentally, it is the second S"*" which belongs to the ground state rotational 
band and our result underestimates its energy by 2.0 MeV. The energy at 
which we obtain a lO"*" state is probably also too low. In the cranking cal- 
culation, the pairing energy becomes very small at spin 8; triaxiality effects 
become also important around this spin. 

In the projected spectra, the energy of the first 2+ state is systematically over- 
estimated, the energy obtained with a volume dependent pairing being how- 
ever rather close to the experimental value. Two restrictions imposed in the 
present calculation may be at the origin of this discrepancy. First, we have not 
included triaxial deformations, which are shown by the cranking calculation 
to be more important at high spins than at low spins. Second, the variation 
after projection on angular momentum is limited to the quadrupole moment; 
the mean-field equations are optimized for the description of the ground state 
energy and not for excited states. The first limitation could be removed with 
minor modifications of the formulae presented in section 2. However, the com- 
puting time would be largely increased, the number of Euler angles for the 
angular integration being at least 50 times larger. A rather simple way to 
improve the variational character of the calculation would be to project for 
each spin wave functions generated by cranking calculations. It is indeed well 
known [8] that cranking is a first order approximation of a variation after 
projection on angular momentum. To generalize our method in this direction 
requires mainly to consider Bogoliubov instead of BCS transformations, a fact 
which would not increase too much the computing time. Work along this line 
is in progress. 

One of the main interests of a restoration of rotational symmetry is the possi- 
bility to calculate transition probabilities without the approximation involved 
in a cranking calculation. On figure 6, are compared to the experimental 
data [33] the transition probabilities along the yrast line obtained in the GCM 
calculation and by considering only the minima of projected energy curves. 
The value of the 2+ to 0'^ B(E2) is nearly independent of the nuclear in- 
teraction and results are only shown for the SLy4 interaction (case b). The 
transition probability between the configurations minimizing the projected en- 
ergy curves is very close to the experimental value. The configuration mixing 
causes a spreading of the collective wave function on the quadrupole moment 
and decreases slightly the value of the B(E2). This effect is similar to the effect 
of quadrupole vibrations that is sometimes included phenomenologically [11] 
in the determination of transition probabilities from intrinsic wave functions. 
Since for spin different from 0, the wave functions do not have components at 
low quadrupole moment, the configuration mixing does not affect significantly 
the transition probabilities. The agreement between both calculations and the 
experimental data is excellent in these cases. 

In our first studies of quadrupole collective dynamics [34], we have introduced 
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an approximate projection on 0^ and 2"*" states from intrinsic states. This 
projection is based on the mixing of the configurations corresponding to all 
the possible labellings of the principal axes of inertia. It is supposed to be 
valid for small deformations. The present calculation gives the opportunity 
to test the validity of this approximate projection. On figure 7 are plotted 
the energy curves corresponding to the mean-field calculation, projection on 
particle number only, projection on particle number together with mixing 
of the three possible orientations of the principal axes and triple projection. 
The effect of particle number projection is dominant for small deformations, 
corresponding to a /3 value up to 0.1. The restricted projection on angular 
momentum gives a fair approximation for slightly larger deformations, up to (3 
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Fig. 6. B{E2) transition probabilities (in e^ fm^ for ^^Mg. Transition probabilities 
between the configurations corresponding to the minimum of the projected energy 
curves of figure 3 (left side) and between the yrast collective states obtained in 
the GCM calculation (right side). In the central part are shown the experimental 
values 
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around 0.2. For still larger quadrupole moments, the full projection becomes 
necessary. On the right part of the figure are compared the configuration 
mixing spectra obtained with the exact and the approximate projections. The 
differences in total energy are due to the failure of the approximate projection 
above quadrupole moments of 0.5 b. However, the relative position of the first 
0^ and 2+ is rather satisfactorily estimated. 

The restricted projection is of limited interest for a well deformed nucleus like 
^^Mg. However for heavy spherical nuclei, its validity up to a /3 value around 
0.2 makes it a cheap alternative to study quadrupole dynamics. 



4 Conclusion 



In this paper, we have presented and tested a method to introduce correlations 
beyond mean-field on HF+BCS wave functions. The formulae given in the first 



10 



> 

CD 



LJJ 








30 



20 



> 



LU 



10 







_io^: 


- 


.6^ ; 


- 


.6^ 


- 


4.+_8^ - 


- 


.0^ 


o:-2* 


2+ 
.4^ 


- 


_6^ 


- 


.2^ 


2! , 
"0^ 


■0^ 




"4^ ; 


-2+ 


-2+ : 


-0^ 


■0^ 


- 



-303-303 
qo(b) 



Fig. 7. On the left part of the figure are shown the energy curves corresponding 
to the mean-field calculation (dashed Une), projection on particle number only (up 
triangles), projection on particle number together with mixing of the three possible 
orientations of the principal axes (down triangles) and triple projection (dots). On 
the right are compared the configuration mixing energies obtained with the full (on 
the left) and the restricted (on the right) projections. 
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part of section 2 are written in a form appropriate for the relaxation of the 
restrictions imposed in this study. The extensive tests performed on the ^^Mg 
nucleus show that the method works with reasonable computing time. 

A first and natural generalization from BCS to full Bogoliubov-Valatin trans- 
formations is under progress. It will allow a better treatment for pairing corre- 
lations. If no significant improvements for even-even nuclei close to the stability 
line is to be expected, HFB is essential to treat correctly nuclei near the drip 
lines (see ref. [31,35]). 

As the projection onto angular momentum implies that one cannot anymore 
make use of the signature symmetry, the generalization to many-body wave 
functions breaking time reversal invariance is a natural next step. That is 
a necessary step towards a description of odd nuclei. It will also enable to 
project wave functions generated for each spin by cranking calculations. As 
it has already been shown theoretically [7,8], the use of cranking wave func- 
tions is the first order of a variation after projection on angular momentum. 
Numerical applications [15,25] have confirmed that the projection of cranking 
wave functions improves the energy obtained for each angular momentum and 
compresses the spectra. Such an effect would correct the too spread spectra 
obtained in the present study. 

Another important question for which we have preliminary answers concerns 
the use of effective interactions adjusted for mean-field calculations in a model 
where are incorporated correlations beyond HF+BCS. The Skyrme parametriza- 
tions that we have tested here behave in a similar way and are reasonable 
starting points to study the effects of correlations. We have started a study of 
several Mg isotopes and of a few neighbouring nuclei is underway. This should 
enable to better determine the properties of effective interactions on which 
depend the spectra of nuclei far from stability. 
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